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I.   INTRODUCTION 

Control  system  complexity  and  demands  on  reliability  are 
rising  with  ever  increasing  automation.  More  and  more  control  system 
engineers  are  looking  to  solid-state  digital  techniques  for  new 
answers  to  the  old  problems  of  reliability,  complexity,  and 
economy.  Digital  techniques  provide  many  advantages  such  as  accuracy 
and  computational  speed  as  well  as  flexibility  and  versatility, 
which  analog  techniques  could  not  provide.  By  the  use  of  digital 
techniques  in  control  systems,  sensitive,  and  efficient  control 
devices  can  easily  be  built.  Another  significant  advantage  of 
digital  control  systems  is  that  the  utilizing  of  digital  components 
in  a  control  system  allows  time  sharing  of  important  parts 
of  the  system  which  results  in  economy  in  the  use  of  equipment. 
The  digital  elements  also  require  much  less  power  than  the  analog 
elements.  The  data  signals  in  digital  control  systems  can  also  be 
transmitted  and  received  in  pulse-code  form  almost  error-free, 
even  in  transmission  through  a  noisy  media.  Furthermore,  a 
digital  computer  in  control  gives  another  important  adventage,  it 
makes  feasible  system  compensation  by  nonlinear  programming  and 
adaptive  or  self -optimizing  control. 

The  digital  servo  system  is  becoming  widely  accepted  in 
machine  tool  control,  radar  systems,  navigational  system,  and 
many  other  fields.  A  digital  servomechanism  uses  numerical 
information  to  provide  a  simple  but  highly  accurate  method  of 
positioning  control.  Figure  1  shows  a  block  diagram  of  a  typical 
digital  servo  system.   The  major  components  of  the  basic  digital 
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servo  system  are  : 

1.  Controlled  system 

2.  Control  elements 

3.  Analog-to-Digital  converter 
**.  Digital-to-Analog  converter 

5.  Digital  data  processing  equipment  or  computer 
The  controlled  system  is  normally  a  linear  servomotor  but  step 
motors  have  been  used  for  quite  some  time  in  both  closed  loop  and 
open  loop  digital  positioning  systems.  The  digital  computer  is  used 
as  the  central  control  unit  for  data  processing,  system  compensa- 
tion, and  decision  making  in  some  cases.  In  a  simple  digital  servo 
system,  a  comparator  such  as  an  up-dovm  counter  can  be  used  . 
The  digital-  to-analog  and  analog-to-digital  converters  are  used  to 
interconnect  the  digital  and  analog  quantities  for  controlled  and 
control  elements.  One  familier  type  of  analog -to -digital  converter- 
encoder  is  the  shaft-position  encoder  disk  employing  a  gray  or 
binary  code.  This  digitizer  is  used  to  translate  the  shaft 
position  into  digital  notation  by  making  use  of  a  mechanically 
rotating  shaft  with  coding  disk  centered  on  the  shaft.  Figure  2 
illustrates  a  typical  coding  disk  for  a   shaft    digitizer.   The 

brushes  contact  the  disk  along  the  reading  line   at   the    point 
shown  with  a  single  brush  positioned   in   the   center   of   each 

concentric  band.  These  concentric  rings  represent  2,  2,  2,. ..,2, 

and   each  ..  ring.'  is  .constructed   of   several  segments  made   of 

either  conducting  material  (the  darkened  areas)   or    insulating 

matorial  (  the  unshaded  areas).  A  signal  is  connected   to    the 
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conductor,  and  if  a  given  brush  makes  contact  with  a  segment  of 
conducting  material,  a  »  1  "  signal  will  result,  but  if  the  brush 
is  over  the  insulating  area,  the  output  from  the  brush  will  be 
a  "  0  ".  The  resolution  of  this  encoder  is  determined  by  the 
number  of  channels.  For  a  ^-channel  coding  disk,  a  resolution  of 
one  part  in  sixteen  is  achieved.  Two  types  of  encoder  disk,  the 
gray  and  binary  code  disks  are  shown  in  Figure  3.  The  coded  disk 
may  be  made  of  transparent  and  opaque  materials  arranged  in  coded 
patterns.  The  photocells  and  sources  are  used  to  read  the  angular 
position  in  a  digital  code. 

In  operation,  the  analog-to-digital  converter    converts    the 
analog  output  of  the  controlled  system  into  a  numerical  quantity, 
which  is  compared  with  a  numerical  order,  and  the  difference  is 
converted  by  the  digital-to-analog  converter   into  an  analog  signal 
capable  of  driving  a  servomotor  until  the  difference  is  eliminated. 

A.   OBJECTIVES 

The  methods  for  analysis  and  synthesis  of  the  analog  system 
such  as  Bode  plots,  Root  Locus,  Parameter  plane  plots,  etc.  have 
been  established  for  a  long  time.  These  methods  work  very  well, 
however,  for  the  discrete  or  digital  system,  there  is  no  such 
method  that  is  as  good  as  those  for  the  analog  system.  Even 
the  well  known  z-transform  method  can  become  very  complicated  if 
the  order  of  the  systems  analyzed  are  higher  than  third  or  fourth 
order. 

The  objective  of  this   paper   was   to   invectigate  a  typical 
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digital  positioning  servo  system  in  order  to  obtain  some  good 
methods  for  analysis  and  synthesis.  In  the  proposed  system,  light 
travels  from  a  light  source  through  equally  spaced  openings  in  the 
code  disk  to  activate  the  sensor.  The  sensor,  in  turn,  generates 
feedback  pulses  at  equilibrium  angular  position.  The  sampling 
period  of  this  system  is  not  constant,  depending  on  the  velocity 
of  the  motor,  the  z-transform  method  of  analysis  and  design  cannot 
be  used  here.  Another  objective  was  to  find  a  suitable  mathematic 
representation  of  such  a  system. 

All  of  the  studies  had  been  done  by  simulation  using  Digital 
Simulation  language  DSL/360  on  the  IBM  360  computer  system,  this 
language  made  the  simulation  quite  easy,  however,  care  must  be 
taken  in  order  to  obtain  accurate  integration  in  the  system  with 
high  velocity. 

In  Chapter  II,  the  system  identification  and  method  of 
simulation  are  given.  Chapter  III  presents  the  discussion  of  the 
results  obtained  from  the  simulation  of  second  and  third  order 
servo  systems.  This  includes  the  study  of  the  stability  limit, 
percent  overshoot,  and  the  effects"  of  variation  of  the 
servoamplifier  gain.  This  chapter  also  shows  the  results  obtained 
from  investigating  the  behavior  of  the  system  with  ramp  input. 
The  Bode  plot  and  parameter  plane  plot  were  used  to  compensate 
the  analog  system  in  Chapter  IV.  The  results  were  then  used  in 
simulating  the  digital  system,  since  there  was  no  direct  compensation 
method. 

The  results  were  obtained  both  in  graphs  and  printed  output, 
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however,  only  graphs  are  included.  Listing  of  all  computer  programs 
are  presented  after  the  conclusion  and  recommendation  in  Chapter  V. 
In  each  study,  the  digital  system  with  several  numbers  of 
opening  slots  and  the  continuous  counterpart  were  simulated.  The 
studies  of  the  digital  systems  were  based  on  the  results  of  studies 
of  the  continuous  system. 
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II.   THE  DIGITAL  POSITIONING  SERVO 
A.  PHYSICAL  PROPERTIES 

The  digital  positioning  servo  system  under  consideration  consists 
of  a  digital  comparator,  a  digital-to-analog  converter,  a  controlled 
system,  and  an  analog -to-digital  converter.  The  analog-to-digital 
converter  consists  of  a  low  inertia  encoder  disk  mounted  tightly 
on  the  servomotor  shaft,  two  photoelectric  sensor  units,  and  a 
direction-sensing  logic  circuit.  The  encoder  disk  has  uniformly 
spaced  slots  around  its  circumference.  Each  photoelectric  sensor 
unit  consists  of  a  light  source  and  a  photodiode.  Figure  it- 
illustrates  the  encoder  disk  and  the  sensor  units  and  Figure  5 
illustrates  the  direction-sensing  circuit. 

Each  time  an  opening  passes  through  the  beam  from  the  light 
source,  a  pulse  is  generated  by  the  sensing  unit,  and  it  will 
eventually  be  fed  to  the  up-down  counter  which  acts  as  the 
digital  comparator.  The  second  sensor  unit  is  for  the  purpose  of 
direction  sensing.  When  the  servomotor  turns  in  forward  direction, 
the  feedback  pulse  becomes  a  down-enable  signal,  when  the  servomotor 
turns  backward,  the  feedback  pulse  becomes  an  up-enable  signal. 
The  two  sensors  are  adjusted  so  that  their  respective  outputs,  X 
and  Y,  overlap  as  shown  in  Figure  6.  Figure  6-A  and  6-B  show 
the  inputs  and  outputs  of  the  direction-sensing  circuit  when  the 
servomotor  turns  forward  and  backward,  respectively.  Figure  7  shows 
a  typical  digital-to-analog  converter. 

For  an  input  as  a  stop  function,  an  appropriate  binary  number 
representing  the  size  of  the  step  function  is  fed  into  the  up-down 
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counter.  The  output  of  the  up-down  counter  is,  in  turn,  converted 
into  an  analog  signal  and  is  sent  to  the  servoamplifier  to 
drive  the  servomotor.  The  sensing  device  will  measure  the  shaft 
position,  encode,  and  send  this  information  to  the  up-down  counter. 
The  positioning  process  will  be  finished  when  the  binary  number 
in  the  up-down  counter  becomes  zero. 

3.  COMPUTER  SLXUUTION 

The  DSL/360  language  was  specially  designed  for  use  in 
simulation.  DSL/3^0  provides  many  advantages  over  other  languages 
and  subroutines  in  existence,  especially  in  simplicity  and  time 
saving.  DSL/360  not  only  provides  these  advantages  in  simulating 
the  linear  physical  system,  but  in  simulation  of  the  non-linear 
system  also.  DSL/3^0  provides  a  basic  set  of  function  blocks 
from  which  a  physical  system  may  be  modeled  such  as  integrators, 
quantizer,  pulse  generators,  function  generators,  and  so  on.  In 
addition,  FORTRAN  Library  functions  from  the  Scientific  Subroutine 
Package  may  be  utilized.  The  user  may  also  provide  his  own 
subroutines  if  necessary. 

In  programming,  first,  each  part  of  the  digital  positioning 
servo  system  was  separated  into  single  function  blocks  as  shown 
in  Figure  8.  Next,  each  single  function  block  was  replaced  by  a 
related  DSL  function  or  standard  FORTRAN  Library  function  as 
appropriate.  Finally,  all  the  DSL  and  FORTRAN  Library  function  were 
put  together  in  the  DERIVATIVE  REGION  of  the  program,  One  advantage 
obtained  in  using  the   DERIVATIVE  REGION  is  that  the  functions  are 
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not  required  to  be  in  any  specific  order. 

In  the  early  parts  of  simulation.  The  servomotor  was  assumed 
to  be  a  second  order  motor  with  the  Transfer  function 

Kv 

Gjs)  =  . 

m     s(s  +  1) 

Later,  the  third  order  servomotor  with  the  transfer  function 

Kv 


G  (s)  = 

m     s(s  +  l)(0.2s  +  1) 

was  used.  The  FORTRAN  Library  function 

Y  =  AINT(X) 
was  used  to  modelled  the  quantizer.   This  function  was  also   used 
to  quantize  the  ramp  input  function.   The  gain  of  the   digital- 
to-analog  converter  was  designed  to  be  equal  to  l/K  to  make  the 
system  a  unity  feedback  system. 

The  fourth  order  Runge-Kutta  with  fixed  integration  interval 
(INTEG  RK3FX)  Integration  routine  was  used  in  all  simulations  to 
assure  the  accuracy  of  the  results.  It  is  very  important  to 
choose  a  correct  integration  interval  (DELT),  since  too  large  DELT 
will  cause  inaccurate  integration  and  too  small  DELT  will  result 
in  excessive  computing  time.  The  DELT  used  in  each  simulation 
is  approximately  five  to  ten  percent  of  the  time  interval  between 
two  consecutive  pulses  generated  by  the  digitizer  at  its  maximum 
angular  velocity  in  that  specific  simulation. 

In  multiple  runs,  the  TERMINAL  region  was   used   to  terminate 
the  previous  run,  to  provide  the  new  sets  of  parameter  values, 
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and  to  terminate  the  program  after  the  output  data . of  the  last 
run  had  been  plotted.  The  FORTRAN  subroutine  DRWG  was  used  in 
the  SAMPLE  region  to  collect  points  for  each  sample  interval 
(DELS)   to  be  plotted  by  the  CALCOMP  XY  Plotter. 

In  most  cases.  The  digital  system  with  10,  V>t  90,  and  180 
openings  and  the  continuous  system  were  simulated  and  the  results 
were  compared.  A  system  with  one  opening  was  used  in  stability 
limit  studies  since  it  was  the  worst  case.  To  compare  the  outputs 
of  the  digital  system  with  output  of  the  continuous  system,  the 
inputs  of  all  systems  must  be  equal.  If  a  step  input  to  the 
continuous  system  is  1.2566^,  inputs  to  the  digital  system  with 
10,  ^5,  90,  and  180  openings  must  be  2,  9t-18t  and  36  respective- 
ly. 

These  input  values  can  be  readily  calculated  by  the   method 

shown  in  example  below, 
EXAMPLE 

To  calculate  comparable  step  inputs  to  the  continuous  and 
digital  system  with  10  openings. 
From  Figure  8, 

in( continuous)  " 


Vin( digital) 

R 
n 

-  IGBIT 

K 
c 

2 

(R     -  IGBIT) 

OP 
for  zero  initial  conditions, 

in (continuous) 
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2     R 
Vin(digital)   =  "5^~ 

to     have     equal     inputs     to     servomotors,     requires 

2     R 

r  =  — a 

OP 
by  arbitarily  choosing     R  e  2 

R  =  1.25664 
The  values  of  inputs  to  digital  system  with   other   numbers  of 
opening  were  also  based  on  this  selection. 
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III.   DISCU33I0N 

The  first  part  of  this  chapter  contains  a  discussion  of  the 

response  of  the  second  order  servo  system  to  a  step  function 

input.  The  third  order  system  with  both  step  and  ramp  function 
input  is  discussed  in  the  later  part  of  the  chapter. 

A.  STEP  RESPONSE  OF  THE  SECOND  ORDER  SYSTEM 

The  main  concern  in  the  response  of  the  second  order  servo 
system  was  its  stability  limit.  The  forward  gain  of  the  system 
had  been  increased  until  the  stability  limit  was  obtained. 

The  model  of  the  second  order  continuous  system  used  in  this 
simulation  is  always  stable  regardless  of  the  size  of  the  gain. 
The  digital  system  was  also  stable  at  all  value  of  gain  regardless 

of   the   number   of   openings   in  the  position  encoder.  The  gain 

7 

of  the  system  was  increased  to  10   to  make  sure  that  the  system 

would  always  be  stable.  Figure  9i  10 1  H.  12,  13 1  and  14  show  the 
step  response  and  plots  of  the  binary  numbers  in  the  up-down 
counter  versus  time  for  the  digital  system  with  number  of  openings 
equal  to  1,  10,  ^5t  90,  and  180,  and  for  the  continuous  system,  res- 
pectively. Figure  15,  16,  17,  18,  19,  and  20  show  their  respective 
phase  plane  plots.  The  plots  of  the  output  of  the  up-down  counter 
were  not  scaled  correctly,  however,  they  could  be  read  easily. 
Each  quantizied  step  of  the  N,  plot  represents  a  one. 

The  difference  that  can  readily  be  seen  in  the  responses  of 
the  digital  and  continuous  systems  is  the  magnitude  of  both  the 
angular  displacement  and  velocity.  The  difference  is  very  pronounced 
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when  the  number  of  openings  is  small.  The  digital  system  with 
fewer  openings  has  longer  time  delay  between  two  consecutive 
feedback  pulses.  This  time  delay  has  allowed  the  motor  to  move 
farther  with  higher  velocity,  similar  to  the  response  of  the 
system  without  a  feedback  loop,  before  the  excitation  signal  is 
reduced  by  the  next  feedback  pulse.  The  system  with  a  large 
number  of  openings  has  a  very  high  sampling  rate,  therefore  this 
system  behavior  is  much  closer  to  that  of  the  continuous  system, 
a  code  disk  with  20000  openings  has  been  developed.  This  disk 
gives  the  system  almost  the  same  response  as  that  of  a  continuous 
system. 

The  phase  plane  plots  of  the  digital  system  are  not  as 
smooth  as  those  of  the  continuous  system.  They  actually  consist 
of  many  short  sections.  Each  portion  represents  the  delay 
interval  between  two  consecutive  pulses.  This  can  easily  be  observed 
in  the  phase  plane  plot  of  the  system  with  lower  number  of 
openings  such  as  shown  in  Figure  16. 

Figure  21  contains  the  comparison  of  the  angular  displacement 
of  the  analog  and  digital  system  with  10  and  V>  openings.  This 
plot  clearly  shows  that  the  natural  frequency  of  this  second 
order  digital  system  is  smaller  than  that  of  the  analog  system. 
It  indicates  that  the  root  locus  of  the  digital  system  would 
cross  the  CO -axis  closer  to  the  origin  than  the  root  locus  of 
the  continuous  system.  This  distance  becomes  smaller  if  the  number 
of  openings  is  decreased. 
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The  plots  of  the  step  response  of  the  digital  system  also 
consist  of  many  short  sections.  It  is  very  easy  to  observe  this 
phenomenon  in  the  response  of  the  system  with  very  low  gain 
such  as  those  shown  in  Figure  22  and  23.  Figure  2k  demonstrates 
how  this  phenomenon  occurs.  At  the  beginning  of  the  simulation, 
the  system  behaves  like  a  system  without  a  feedback  loop.  The 
rate  of  increase  of  velocity  is  very  high  at  first,  however,  it 
becomes  lower  with  time.  The  displacement  plot  also  becomes  more 
linear  with  time.  When  the  displacement  is  large  enough,  a  feedback 
pulse  is  generated,  the  output  of  the  up-down  counter  is  changed 
by  one,  the  velocity  and  the  rate  of  change  of  the  displacement 
are  changed  to  follow  the  new  input  voltage  to  the  servoamplifier. 
This  process  occurs  repeatly  until  the  output  of  the  up-down 
counter  is  eliminated,  and  the  system  eventually  comes  to  rest. 

Figure  22  and  23  also  illustrate  that  the  digital  system  behaves 
like  the  continuous  system  in  another  aspect.  When  the  gain  is 
reduced,  the  peak  overshoot  decreases.  This  peak  overshoot  can  be 
eliminated  by  making  the  gain  small  enough.  It  is  also  observed 
that,  for  all  cases  studied,  the  steady  state  displacement  of  the 
digital  system  is  always  larger  than  that  of  the  continuous 
system  with  comparable  size  of  step  input.  The  difference,  however, 
has  never  exceeded  the  angular  displacement  between  two  adjacent 
openings.  The  system  with  more  openings  in  the  encoder  disk  gives 
smaller  error  than  the  system  with  fewer  openings.  The  maximum 
error  for  the  system  with  10,  45,  90,  and  180  openings  are 
less  than  0.6283,  0.1396,  O.O698,  and  0,0349  radians  respectively. 
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B.  THIRD  ORDER  SYSTEM 

A  third  order  system  was  studied  in  greater  detail  than 
the  second  order  system.  In  addition  to  what  had  been  studied 
from  the  second  order  system  characteristics,  the  variation  of 
the  percent  overshoot  was  investigated  in  this  part. 

1.  Stability  limit 

A  third  order  continuous  system  with  the  transfer  function 
shown  in  Chapter  II  reaches  the  stability  limit  when  the 
servoamplifier  gain,  K^  is  equal  to  6.  The  digital  systems, 
however,  require  larger  gains  to  reach  their  stability  limits. 
Figure  25,  26,  27,  28,  and  29  illustrate  the  phase  plane  plots  of  the 
third  order  digital  system  with  the  gain  equal  to  6  and  the 
number  of  openings  equal  to  1,  10,  ^5,  90,  and  180  respectively. 
Figure  30  shows  the  phase  plane  plot  of  the  continuous  system 
with  the  same  gain.  The  phase  plane  plots  of  the  third  order 
digital  system  converge  into  a  limit  cycle  and  can  clearly  be 
observed  in  the  case  the  system  has  1  or  10  openings.  By 
inspecting  the  characteristic  of  the  phase  plane  plots,  especially 
in  the  case  where  the  number  of  openings  is  small,  it  is 
obvious  that  the  third  order  digital  system  requires  higher  gain 
to  reach  the  stability  limit  than  the  continuous  system.  The 
time  delay  also  causes  the  displacement,  velocity,  and  limit 
cycle  size  of  the  digital  system  to  be  larger  than  those  of  the 
continuous  system. 

Since  thero  had  not  been  any  known  method  that  can  bo 
used  to  find  the  gain  at  stability  limit  of  this  digital  servo 
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system,  a  trial,  and  error  method  was  used  to  find  this  gain,  the 
gain  of  the  third  order  digital  system  with  10  openings  at  its 
stability  limit  was  found  to  be  approximately  6.25.  This  gain 
would  be  reduced  toward  6  if  the  number  of  openings  is  increased. 
The  size  of  the  limit  cycle  is  also  decreased  with  the  increase 
of  the  number  of  openings. 

The  results  from  simulating  the  third  order  digital  servo 
system  with  number  of  opening  equal  to  10  as  illustrated  in 
Figure  31 »  26,  32,  and  33  also  show  that  the  size  of  the  limit  cycle 
increases  with  the  gain.  The  increase  of  the  gain  also  causes 
the  system  to  reach  this  limit  cycle  faster.  Upon  making  the 
gain  greater  than  the  stability  limit  gain,  The  phase  plane  plot 
diverges  for  a  period  of  time  before  reaching  another  limit  cycle 
as  shown  in  Figure  3^.  For  higher  gain,  The  phase  plane  plot 
diverges  into  a  limit  cycle  and  stays  in  the  limit  cycle  for  a 
period  of  time  then  diverges  to  infinity  as  shown  in  Figure  35. 
With  the  gain  a  little  higher  as  shown  in  Figure  36,  the  phase 
plane  plot  diverges  to  infinity  but  has  a  tendency  to  go  into 
a  limit  cycle  at  first.  When  the  gain  is  much  higher,  the 
system  just  becomes  unstable  as  shown  in  Figure  37. 
2.  Percent  overshoot 

For  a  continuous  and  stable  servo  system,  the  overshoots 
of  the  step  response  always  decrease  with  time,  and  the  plot  of 
the  log  of  the  percent  overshoot  versus  time  is  a  straight 
line.  The  step  response  and  the  plot  of  log  of  percent  overshoot 
of  the  continuous  system  are  shown  in  Figure  38  and  ^3  respectively. 
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To  investigate  the  behavior  of  the  overshoot  of  the  digital 
positoning  servo  system,  a  gain  of  2.5  was  used.  With  this  gain, 
£  and  cjn  of  the  continuous  system  are  equal  to  0.166  and  1.54 
respectively.  Figure  39,  40,  WL,  and  42  illustrate  the  step  response 
and  output  of  the  up-down  counter  of  the  third  order  digital  servo 
system  with  180,  90,  45,  and  10  openings  respectively.  The  step  response 
of  the  continuous  system  shows  6  prominent  overshoots,  the  responses 
of  the  digital  system  with  180,  90,  45,  and  10  openings  show  5,:4, 
3t  and  2  overshoots  respectively.  The  time  at  which  each  of  the 
corresponding  peaks  of  all  systems  occur  is  the  same.  Figure  4-3 
illustrates  the  plots  of  the  percent  overshoot  of  the  responses 
of  the  continuous  and  digital  system  with  180,  90,  45,  and  10 
openings  versus  time.  The  numerical  values  of  percent  overshoots 
are  included  in  table  I.  The  inputs  to  all  systems  were  adjusted 
so  that  they  would  be  equal  and  the  percent  overshoot  could  be 
calculated  from  the  same  basis,  the  steady  state  output  of  the 
continuous  system. 

The  overshoots  of  all  digital  systems  are  higher  than  the 
corresponding  peaks  of  the  continuous  system.  The  difference, 
however,  does  not  exceed  the  angle  of  the  two  adjacent  openings. 
The  digital  system  with  smaller  number  of  openings  has  larger 
percent  overshoot  than  the  system  with  larger  number  of  openings. 
The  steady  state  response  of  the  digital  system  is  also  larger 
than  that  of  the  continuous  system.  The  difference  is  also  less 
than  the  angular  displacement  between  two  adjacent  openings. 

At  the  beginning,  the  plot  of  the  percent  overshoot  of  the 
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digital  system  tends  to  be  parallel  to  that  of  the  continuous 
system.  The  plot  of  the  percent  overshoot  then  diverges  to  the 
smallest  peak  which  is  slightly  larger  than  the  angular  displacement 
between  two  adjacent  openings.  Since  the  system  with  fewer 
openings  has  larger  displacement  between  the  openings,  its  response 
requires  less  time  to  reach  the  last  overshoot.  The  number  of  its 
overshoots  is  also  less  than  that  of  the  system  with  more  openings. 

The  responses  of  all  the  digital  systems  that  have  less 
than  10  openings  were  also  studied.  All  of  the  responses  have 
two  overshoots  and  the  characteristics  are  similar  to  those  discussed 
above.  However,  only  the  response  of  the  system  with  one  opening, 
which  is  the  extreme  case,  is  illustrated  in  Figured. 
3.  Response  to  ramp  input 

The  ramp  input  function,  Bt,  was  also  applied  to  the  third 
order  digital  positioning  servo  system  in  order  to  investigate  its 
stability.  The  servoamplifior  gain  used  in  this  study  was  2,5 
and  the  slope  of  the  ramp  was  varied  from  0.5  to  .5.  .  Practically, 
a  method  of  applying  a  ramp  input  to  the  digital  servo  system 
is  to  apply  a  series  of  impulses  into  the  up-down  counter.  The 
slope  of  the  ramp  can  be  controlled  by  varying  the  period  of  the 
impulse.  A  digitized  ramp  function  was  used  in  this  simulation 
instead  of  the  impulses.  The  study  of  the  response  of  the  digital 
servo  system  with  ramp  input  was  done  by  simulating  the  system 
with  10  openings  only.  It  had  been  shown  from  the  previous  studies 
that  the  digital  system  with  more  openings  behaves  closer  to  the 
continuous  system  than  the  system  with  fewer  openings. 
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Figure  45  and  46  illustrate  the  response  of  the  digital  system 
with  10  openings  and  the  slope  of  the  ramp  input  equal  to  0.5, 
1.0,  1.5,  2.0,  2.5,  3.0,  4.0,  and  5.0.  When  the  slope  is  small, 
for  example  B  =  0.5,  the  rate  of  the  input  pulse  applied  to  the 
up-down  counter  is  slow  and  causes  the  velocity  of  the  shaft  and 
the     feedback     sampling     rate     to     be     too     slow     and     nonuniform.  The 

nonuniform  and  slow  rate  of  output  of  the  up-down  counter,  in 
turn,  cause  the  shaft  to  turn  abruptly  and  the  system  becomes 
unstable.  The  longer  time  delay  is  also  an  important  cause  for 
the  instability  of  the  system.  The  plot  of  velocity  and  output 
of  the  up-down  counter  for  the  digital  system  with  the  slope  of 
the  ramp  input  to  the  system  equal  to  0.5,  i.0,  1.5,  2.0,  2.5, 
3.0,  4.0,  and  5.0  are  shown  in  Figure  47„  48,  49,  5Q,  51,  52,,  53, 
and  54  respectively.  When  the  slope  of  the  ramp  input  is  increased, 
the  input  pulses  are  fed  to  the  up-down  counter  at  higher  rate, 
the  servo  shaft  turns  faster,  the  feedback  sampling  rate  becomes 
higher,  the  output  of  the  up-down  counter  becomes  more  linear  and 
the  velocity  of  the  shaft  becomes  more  stable.  Unlike  the  continuous 
system,  the  velocity  of  the  digital  system  will  not  be  constant 
no  matter  how  large  a  slope  is  used  for  the  ramp  input.  It  rather 
varies  within  a  small  limit.  This  limit  cycle  is  caused  by  the 
time  delay  between  pulses.  The  limit  cycle  will  be  smaller  if 
the     time     delay     is     shorter. 

Another  interesting  result  observed  is  that  the  rate  of 
decrease  of  the  magnitude  of  the  limit  cycle  with  the  increase 
of     the     slope     of     the     ramp     is     nonuniform.        At     first     the  magnitude 
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of  the  limit  cycle  becomes  very  small  when  the  slope  is  increased 
to  2.5.  However,  when  the  slope  is  increased  to  3  and  4-,  the 
magnitudes  of  the  limit  cycles  are  some  what  larger.  The  limit 
cycle  is  reduced  to  a  smaller  size  again  when  the  slope  is 
increased  to  5. 
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IV.  SYSTEM  COMPENSATION 

Since  there  had  not  been  any  known  method  that  can  be  used 
to  compensate  this  digital  positioning  servo  system  directly,  the 
author  chose  to  compensate  the  comparable  continuous  system  with 
known  methods  and  simulated  the  compensated  system  as  both  a 
continuous  and  a  digital  system.  The  results  from  the  simulation 
were  then  compared  and  studied.  The  continuous  and  the  digital 
system  with  10,  k$t  and  180  openings  were  simulated  and  their 
step  responses  were  plotted  on  the  same  graph  to  make  the 
comparison  more  obvious.  The  simulation  of  the  digital  system  with 
90  openings  was  omitted  to  prevent  the  graph  from  being  overcrowded. 
The  system  compensated  was  the  third  order  servo  with  servoamplifier 
gain  equal  to  6. 

A.  LAG  COMPENSATION 

Since  the  programs  and  computer  facility  were  available,  the 
parameter  plane  plot  method  was  used  in  lag-compensating  of  the 
system.  This  method  is  very  efficient  and  flexible  and  the  computer 
make  it  very  easy  to  obtain  the  plot.  The  author  selected  the 
PARAl'  M  program  to  use  in  this  compensation.  The  complete  details 
of  this  program  are  shown  in  the  appendix. 

In  order  to  obtain  the  parameter  plane  curve,  the  characteristic 
equition  of  the  system  and  other  requirements  such  as  the  value 
of  C  »  CJn  »  Co  ,  or  <3  needed  to  be  plotted  and  the  scales  of  the 
plot  must  be  provided  as  the  input  data  to  the  PARAK  M  program. 
Figure  55  illustrates  the   block    diagram   of   the   third   order 
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continuous   system   with   the  lag  compensation   section.   The  lag 
compensation  block  has  the  transfer  function 

(T„S  +  1) 


G  (s)  = 


(Tps  ♦ 1) 


The  transfer  function  of  the  third  order  system   with   the   lag 
compensation  becomes 

(t>  +  i)        6 


G  (s)  = 


m 


(TpS  +  1)  s(s  +  l)(0.2s  +  i) 
or 

g(s)  =El!_til 2° 

z  (s  +  p)  s(s  +  l)(s  +  5) 

The     characteristic     equation     of     this     system     is 

s^  +  (6  +  p)s3  +  (5  +  6p)s2  +  (5P  +  3o|  )s  +  30p  =  0 
let  oc  =  p 

the     characteristic     equation     becomes 

s^  +  (6  +a)s3  +  (5  +  6oc)s2  +  (5  +  30p)s  +  30oc=  o 
The  parameter  plane  plot  of  this  characteristic  equation  is  shown 
in  Figure  5^.  cc  and  £  are  the  unknown  pair  of  variables  which 
their  values  can  be  directly  selected  from  the  suitable  value  of 
£  and  GJn.  The  values  of  pole  and  zero  can  also  be  readily 
obtained  from  the  values  of  oc  and  0 .  Three  pairs  of  poles  and 
zeroes  were  chosen  and  simulated.  The  step  responses  of  the 
compensated     continuous     and     digital     system     with     10,        ^5,      and       180 
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openings  are  illustrated  in  Figure  57,  58,  and  59.   The  values  of 

pole,   zero,  C,  6Jn ,  and  predicted  and  resulted  values  of  M  ..  t  . 

pt'  p* 

and  t   of  the  compensated  system  are  included  in  table  II. 

The  simulated  results  show  that  the  characteristics  of  the 
digital  system  are  similar  to  those  of  the  continuous  system. 
The  results  also  agree  with  those  studied  in  the  previous  chapter. 
The  magnitude  of  the  peak  overshoot  and  steady  state  response  of 
the  digital  system  are  always  larger  than  those  of  continuous 
system.  All  the  differences  are  also  within  the  limit  of  the 
angular  displacement  between  two  adjacent  openings. 

Care  must  be  taken  in  selecting  the  value  of  the  zero  for 
the  lag  compensation  section.  If  this  zero  was  too  close  to  the 
origin,  the  real  root  would  dominate  the  system  instead  of  the 
selected  pair  of  imaginary  roots.  If  the  zero  was  too  far 
from  the  origin,  the  value  of  £  would  be  too  small  and  cause 
the  system  to  be  very  slow. 

B.  LFAD  COMPENSATION 

The  same  character jstic  equation  was  used  to  obtain  the  parameter 
plane  plot  for  lead  compensation.  However,  the  scales,  the  starting  value 
of  COfl  ,  "the  required  values  of  5  a"d  U>n  curves,  and  etc.  must 
be  res elected  in  order  to  move  the  curves  into  required  region 
for  the  lead  compensation.  The  block  diagram  of  the  system  with 
lead  compensation  section  is  the  same  as  that  of  the  lag 
compensated  system.  Its  parameter  plane  plot  is  shown  in  Figure  60. 

The  pole  and  zero  of  the  lead  compensation  section  were  kept 
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within  one  decade  apart  to  prevent  the  severe  reduction  of  the 
forward  transfer  gain  of  the  section.  This  imposes  another  limitation 
into  the  selection  of  the  value  of  (3.  The  value  of  £  should 
not  be  larger  than  10.  Three  pairs  of  poles  and  zeroes  were 
selected  and  simulated.  Figure  61,  62,  and  63  illustrate  the  step 
response  of  the  compensated  system  in  both  continuous  and  digital 
versions.  Table  IU  consists  of  the  values  of  pole,  zero,  £  ,  cjn, 
predicted  and  resulted  values  of  t  ,   M  . ,  and  t  of  the  compensated 

system.   Since  the  natural   frequency   of   the   lead   compensated 

system  is  larger  than  that  of  the   lag   compensated   system   it 

results  in  smaller  t   and  t  .  The  other  results  of  this   study 

p       s  J 

are  similar  to  those  described  in  the  lag  compensation  study  above. 

C.  LAG -LEAD  COMPENSATION 

The  block  diagram  of  the  third  order  continuous  system  with 
lag  and  lead  compensation  sections  is  shown  in  Figure  6k,  There 
are  four  variables  in  this  system  which  make  it  more  difficult 
to  use  the  parameter  plane  plot  method  in  compensation.  The  Bode 
plot  method  was  used  in  this  compensation.  This  method  is  rather 
simple  but  requires  some  degree  of  neatness  .  Figure  65  a"d  66 
illustrate  the  Bode  plots  of  two  different  lag-lead  compensations.. 
The  limitation  in  selecting  the  values  of  the  zero  of  the  lag 
section  and  the  spread  between  the  pole  and  the  zero  of  the  lead 
section  are  also  imposed  here.  Both  of  the  compensated  systems 
were  also  simulated  and  studied.  The  results  of  the  simulation 
of  the  lag-lead  compensated  systems  are  shown  in  Figure  6?  and  68, 
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Table  IV  consists  of  the  values  of  pole,   zero,  5,  (Jr\ »   predicted 
and  resulted  value   of  t  ,   M  .  ,  and  t  .   The  results  of  the  lag- 

P     pT*  S 

lead  compensation  study  are  also  similar  to  those  of  lag  and  lead 

compensation  studies. 

D.  COMPENSATION  DISCUSSION 

The  results  from  all  compensation  studies  as  well  as  those  from 
the  previous  chapter  show  that  the  steady  state  response  of  the 
digital  system  is  always  within  the  angular  displacement  between 
two  adjacent  openings  larger  than  that  of  the  comparable  analog 
system.  However,  this  difference  is  not  constant,  it  changes  randomly 
within  the  limit.  The  exact  position  of  the  steady  state  response 
of  the  system  can  not  be  determined.  It  seems  that  the  best  method 
to  compensate  the  digital  system  is  to  compensate  its  analog 
counterpart. 

The  maximum  error  of  the  system  under  investication  is  found 
to  be 

1.296X1Q6    seconds 
OP 

In  comparison,   the  analog  system  has  the  accuracy  of  approximately 

2  minutes   of   arc   and   the   digital   system  with  gray  or  binary 

shaft-position  encoder  has  the  accuracy  of 

1.S*X106    seconds 
2n 
where  n  is  the  number  of  the  conducting  channel  (see  ref.  5).   In 
order  to  have  the  same  accuracy  as  of  the  analog   system,   this 
digital  system  must  have  approximately  iOOOO  openings  in  the  encoder 
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disk.  The  digital  system  with  gray  or  binary  code  disk  requires 
13  channels  in  the  code  disk  to  acquire  comparable  accuracy.  It 
is  trivial  that  this  digital  system  has  an  economic  advantage  over 
the  gray  or  binary  code  disk  system.  This  system  requires  only 
two  sensing  units  for  any  number  of  openings.  The  gray  or  binary 
shaft-position  encoder  unit  requires  one  sensing  device  for  each 
channel  in  the  code  disk  and  the  sensing  circuit  has  to  be  more 
complicated. 

Another  adventage  is  that  all  the  studies  show  that  the 
response  of  this  digital  system  is  always  larger  than  the  response 
of  the  continuous  system.  This  observation  could  be  used  to  reduced 
the  steady  state  error  by  making  the  size  of  the  step  input  to 
be  less  than  the  required  input  to  the  system.  The  difference 
should  be  as  close  to  a  half  of  the  maximum  error  as  possible. 
The  difference  of  exactly  a  half  of  the  angular  displacement 
between  two  adjacent  openings  would  reduced  the  steady  state  error 
to  within  one  half  of  the  maximum  error. 
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■  C0NCLU33I0N 

The  characteristics  of  the  digital  positioning  system  under 
investigation  can  be  very  similar  to  the  analog  system  if  the 
feedback  sampling  rate  is  high  enough.  This  system  has  some 
characteristics  which  are  superior  to  the  conventional  analog  and 
sampled-data  system,   especially  in  accuracy  and  economy. 

Compensation  of  this  system  can  be  accomplished  by  using  the 
results  obtained  from  compensating  its  analog  counterpart  by 
conventional  methods.  The  accuracy  can  be  readily  improved  by 
modifying  the  encoder  disk. 

From  the  discussion  above,  it  can  be  seen  that  there  are 
many  areas  left  to  be  studied.  The  author  believes  that 
further  study  of  the  mathematical  analysis  and  of  the  direct  method 
of  compensation  are  necessary. 
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Figure  1.   Block  diagram  for  a  typical  digital  control  system. 
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Figure  2.  A  typical  shaft  position-to-digital  converter 
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Figure  15.  Second  order  digital  servo,  K  =10  ,  R  =1,  0P=1 

phase  plane  plot 
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Figure  16.      Second  order  digital  servo,   Kv=10   ,   ^n=2»   OP=10 

phase  plane  plot 
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Figure  17.  Second  order  digital  servo,  K  =10  ,  R  =9,  0P=45 


phase  plane  plot 
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Figure  18.     Second  order  digital  servo,   K  =1(T  t   R  =18,   OP=90 

phase  plane  plot 
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Figure  19.  Second  order  digital  servo,  Kv=10  ,  ^n=36,  CP=180 

phase  plane  plot 
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phase  plane  plot 
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Figure  22.   Step  response  of  second  order  digital  servo 
Kv=1.5,  1.0,  0.5,  0.36,  0.25,  0.1,  and  0.05, 
0P=10,  and  R  =2 
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Figure  23.  Step  response  of  second  order  digital  servo 
K  =1.5,  1.0,  0.5,  0.36,  0.25,  0.1,  and  0.05, 
0P=45,  and  R  =2 
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Figure  2*f.   Second  order  digital  servo,  K  =0.05,  OP=10,  and  R  =2 
°  °  v  n 

Up-down  counter  output,  velocity,  and  displacement  vs  time 
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Figure  25.   Third  order  digital  servo,  K  =6,  R  =1,  0P=1 

phase  plane  plot 
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Figure  26.  Third  order  digital  servo,  K  =6,  R  =2,  OP=10 

phase  plane  plot 
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Figure  27.  Third  order  digital  servo,  K  =6,  R  =9,  OP=4-5 

phase  plane  plot 
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Figure  28.  Third  order  digital  servo,  K  =6,  R  =18,  OP=90 

v     n 

phase  plane  plot 
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Figure  29.  Third  order  digital  servo,  K  =6,  R  =36,  OP=180 

phase  plane  plot 
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Figure  30.   Third  order  analog  servo,  K  =6,  R=l.  25664- 

phase  plane  plot 
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Figure  31.  Third  order  digital  servo,  Kv=5.7,  R  =2,  OP=10 

phase  plane  plot 
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Figure  32.      Third  order  digital  servo,   K  =6.1,   R  =2,   0P=10 

v  n     * 

phase  plane  plot 
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Figure  33.   Third  order  digital  servo,  Kv=6.25,  R  =2,  OP=10 

phase  plane  plot 
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Figure  3^.   Third  order  digital  servo,  K  =6.5,  R  =2,  0P=10 

phase  plane  plot 
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Figure  35.   Third  order  digital  servo,  Kv=7.1,  Rn=2,  OP=lO 

phase  plane  plot 
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Third  order  digital  servo,   K  =7.2,   R  =2,   OP=10 
phase  plane  plot 
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Figure  37.  Third  order  digital  servo,  K  =7.5,  R  =2,  0P=10 

v      n 


phase  plane  plot 


78 


-  « 


s 

2 


© 

E 
© 
o 
OS 
iH 
D. 
W 


— 1  /     \ 

p               \          1 

• 

tr> 

<*i 

a 
i 

rs/ 

i 

«■- 

•a 

1                                    1 

1                                    1                                   1 

6-00 


12.  ca 


10.  00 


24.00 


30 


Time,   seconds 

Figure  38.   Third  order  analog  servo,  Kv=2.5,  R=l. 25664- 

displacement  versus  time 
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Figure  39.      Third  order  digital  servo,   K  =2.5,   R  =36,   OP=180 

v       n 

displacement  and  up-down  counter  output  vs  time 
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Figure  *+0.  Third  order  digital  servo,  K  =2.5,  ^.=18,  OP=90 
displacement  and  up-down  counter  output  vs  time 
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Figure  M..  Third  order  digital  servo,  K  =2.5,  R  =9,  0P=45 
displacement  and  up-down  counter  output  vs  time 
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Figure  hZ.      Third  order  digital  servo,   K  =2.5,  R  =2,  OP=10 
displacement  and  up-down  counter  output  vs  time 


63 


100- 


o 
o 

« 

o 
> 

o 

c 

© 


c 
&, 


50 


: 


10 


5  - 


1- 


0.5- 


a  OP  =  10 

x  OP  =  k$ 

o  OP  s  90 

+  OP  =  180 

A  ANALOG 


10 
Time, 


15 

seconds 


20 


Figure  kj.  Third  order  servo  system,  K  =2.5 
OP=10,  ^5,  90,  180,  and  analog 
percent  overshoot  versus  time 
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Figure  W*.   Third  order  digital  servo,  K  =2.5,  R  =2,  0P=1 


displacement  versus  time 
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Figure  ^5.  Third  order  digital  servo,  displacement  vs  time 
ramp  input,  B=0.5,  1.0,  1.5, and  2.0,  OP=10,  K  =2.5 
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Figure  ^-6.  Third  order  digital  servo,   displacement  vs  time 
ramp  input,  3=2.5.  3.0,  ^.O.and  5.0,,0P=10,  Ky=2.5 
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Figure  b?.   Third  order  digital  servo,   ramp  input,  B=0.5,   Kv=2.5,   OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  48.  Third  order  digital  servo,  ramp  input,  B=1.0,  Kv=2.5,  0P=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  *J-9.  Third  order  digital  servo,  ramp  input,  3=1.5,  K  =2.5,  OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  50.  Third  order  digital  servo,  ramp  input,  B=2.0,  K  =2.5,  OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  51.    Third  order  digital  servo,   ramp  input, .B=2.5,   Kv=2.5,   OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  52.  Third  order  digital  servo,  ramp  input,  B=3.0,  K  =2.5,  OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  53.  Third  order  digitals ervo,  ramp  input,  B=4-.0,  Kv=2.5,  OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  5^.  Third  order  digital  servo,  ramp  input,  B=5.0,  Kv=2.5,  OP=10 
velocity  and  up-down  counter  output  versus  time 
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Figure  56.  Parameter  plane  plot  for  lag  compensation 
of  the  third  .order  positioning  servo 
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Figure  57.   Lag  compensated  third  order  positioning  servo 
p=0.0!3,   a=0.18,     0P=10,   ^5,    ^80,   and     analog 
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Figure  58.  Lag  compensated  third  order  positioning  servo 
P=0.017,  z=0.19,  0P=10,  ^5,  180,  and  analog 
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Figure  59.  La?  compensated  third  order  positioning  servo 
pcrO.Ol,  £=0.15,  0P=10,  ^5,    130,  and  analog 
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Figure  60.  Parameter  plane  plot  for  lead  compensation 
of  the  third  order  positioning  servo 
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Figure  61.  Lead  compensated  third  order  positioning  servo 
P=i4.5,  z=1.85,  0P=10,  ^5 1  180,  and  analog 
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figure  62.  Lead  compensated  third  order  positioning  servo 
P=15.0,  ^=1.50,  OP=10,  ^5,  180,  and  analog 
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Figure  63.  Lead  compensated  third  order  positioning  servo 
p=l9.0,   2=1.90,  OFrrlO,  ^5,  180,  and  analog 
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Figure  6?.   Lag-lead  compensated  third  order  servo 

p_      =0.25,   z_      =0.5,   p,       ,=20,   z.       ,=2 
*lag  lag       -"    rlead       *     lead 

OP  =  10,    ^5,    180,     and       continuous 


108 


i 

<D 
O 

•5 

a. 
n 


4.23 


6-0: 


12. 


Time,   seconds 


1S-QQ 


2( 


Figure  68.   Lag-lead     compensated     third     order     servo 

Plag=°-2'   ZlaB=°-5'  Pload=12-5.   *lead=1-25 

OP  =  10,   4-5,    180,     and     continuous 
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APPENDIX 

This  appendix  presents  four  examples  of  the  programs  used  to 
simulate  the  system  under  investigation  and  a  program  used  to  plot 
the  parameter  plane  curves. 
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TITLE       SECOND  ORDER  ANALOG  POSITIONING  SERVO 

*  STEP  RESPONSE  AND  PHASE  PLANS  PLOT 

* 

* 

INTEG  RK3FX 

INTGER  NPLOT.NCUR 

CONST  NPLOT=2,   NCUR=2,KV=1.0E7,R=1.2566^ 

INITIAL 

DERIVATIVE 

VT=R-C 

CDD=VT*KC 

CDOT=REALPL(0.0,l.O,CDD) 

C=INTGRL(0.0,CDOT) 
CONTRL  FINTIX=O.Ol,DSLT=0. 00000 1,DELS=0. 000002 
PRINT  0.000l,VTfCDOT,C 
SAMPLE 

CALL  DR>/G(i,lTLT2,C) 

CALL  DR'.VG(2,l,CfCD0T) 
TERMINAL 

CALL  ENDRW(NPLOT) 
END 
STOP 
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TITLE       SECOND  ORDER  DIGITAL  POSITIONING  SERVO 

*  STEP  RESPONSE  AND  OUTPUT  OF  THE  UP-DOWN  COUNTER 

* 

*  A  IS  A  FACTOR  USED  TO  ADJUST  THE  SIZE  OF  THE  CURVE 

*  MULTIPLE  RUNS 
* 

INTEG  RK3FX 

INTGER  NFLOT.NCUR 

CONST  ^LOT=1,NCUR=5,OP=1,RN=1,KV=1.0E7,A=20.0 

INITIAL 

KC=OP/6. 2831853 
DERIVATIVE 

BIT=C*KC 

IGBIT=AINT(BIT) 

NT=RN-IG3IT 

VT=NT/KC 

CDD=VT*KV 

CDOT=RSALPL( 0.0,1.0, CDD) 

C=INTGRL(0.0,CDOT) 

N0=NT*A 
CONTRL  FTNTBfcO  .01 ,  DSLT=0 .  OOOOOl ,  DELS=0. 00002 
PRINT  0.0001,VT,CDOT,C,NT 
SAMPLE 

CALL  DRWG(NPL0T,1,TIMS,N0) 

CALL  DRV/G  ( NPLOT ,  2 ,  THE ,  C ) 
TERMINAL 

IF(NCUR.SQ.l)   GO  TO  6 

NPLOT=NPLCT+l 

NCUR=NCUR-1 

IF(NCUR.EQ.^)   GO  TO  1 

IF(NCUR.E2.3)   GO  TO  2 

IF(NCUK.BQ.2)   GO  TO  3 

IF(NCUR.EQ.l)   GO  TO  4 

1  OP=10.0 
RN=2.0 
A=1.5 
GO  TO  5 

2  OP=45.0 
RN=9.0 
A=0.3 
GO  TO  5 

3  op=90.o 

RN=18.0 
A=0.15 
GO  TO  5 
^  OP=180.0 
RN=36.0 

A=0.075 

5  CALL  RERUN 
RETURN 
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6  CALL  ENDRW(NPLOT) 


end 

STOP 
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TITLE  THIRD  ORDER  DIGITAL  SERVO 

*  VELOCITY  AND  OUTPUT  OF  UP-DOWN  COUNTER 

*  VARIATION  OF  SLOPS  OF  RANP  FUNCTION 

* 

*  MULTIPLE  RUNS 

INTEG  SK3FX 

INTGER  NPLOT.NCUR 

CONST  NPLOT=l,NCUR=4,KV=2. 5, OP=10. 0,3=5.0 

INITIAL 

KC=0P/6. 2831 853 
DERIVATIVE 

RAK=3*RAI?(0.0) 

RN=Arrr(RAM) 

BIT=C*KC 

IG3IT=AINT(BIT) 

NT=RN-IG3IT 

VT=NT/KC 

CD=VT*KV 

CDD=REALPL(0.0,0.2,CD) 

CD0T=REALPL  (  0 . 0 , 1 . 0 ,  C  DD) 

C=INTGRL(0.0,CDOT) 

NO=NT 
CONTRL  FINTT-^O .  0 ,  DELT=0 . 004,  DELS=0 .03 
PRINT  0.2,VT,CD0TtC,NT 
SA1-TLE 

CALL  DRWG ( NPLOT ,  1 , TIME , NO ) 

CALL  DRV/G(NTLOT, 2,Ti:-S, CDOT) 
TERMINAL 

IF(NCUR.EQ.l)   GO  TO  6 

NCUR=NCUR-1 

NPL0T=NPL0T+1 

IF(NCUR.E3.3)   GO  TO  1 

IF(NCUR.EQ.2)   GO  TO  2 

IF(NCUR.EQ.l)   GO  TO  3 

1  3=4.0 
GO  TO  5 

2  3=3.0 
GO  TO  5 

3  B=2.5 

5  CALL  RERUN 
RETURN 

6  CALL  ENDRW(NPLOT) 
END 

STOP 
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TITLE.      THIRD  ORDER  POSITIONING  SERVO 

*  STEP  RESPONSE 

* 

*  MULTIPLE  RUNS 

* 

*  ANALOG  AND  DIGITAL  SYSTEMS 

INTEG  RK3FX 

INTGER  NPLOT.NCUR 

CONST  ^LOT=4,::CUR=l,C=O.OfD=1.0,KV=6.0,RN=2.0,OP=lO.O,Pl=0.2,... 

Z1=0.5»P2=12.5,Z2=1.25 
INITIAL 

KC=OP/6. 283 1853 

TZ1=1.0/Z1 

TPl=1.0/Pl 

TZ2=1.0/Z2 

TP2=1.0/P2 
DERIVATIVS 

BIT=KC*C*D 

IGBIT=AINT(BIT) 

NT=RN-IG3IT-C*C 

VT=NT*(D/KC+C) 

VTN=VT*KV 

CD=LSDLAG(O.OfTZl,TPl,VTN) 

COD=LEDLAG(0.0,TZ2,TF2,CD) 

CDD=REALPL(0.0,0.2,COD) 

CDOT=REALPL(0.0,1.0,CDD) 

C=INTGRL(0.0,CDOT) 
CONTRL  FINTiys=20.0,DSLT=0. 0025, DELS=0. 015 
PRINT  0.2,VTfCD0T,C,NT 
SAMPLE 

CALL  DRV/G(1,NCUR,TIKE,C) 
TERMINAL 

IF(NPLOT.EQ.l)  GO  TO  6 

NPLOT=NPLOT-l 

NCUR=NCUR+1 

IF(KPLOT.E3.3)  GO  TO  2 

IF(NPL0T.EQ.2)  GO  TO  3 

IF(NPLOT.EQ.l)  GO  TO  4 

2  RN=9.0 
OP=45.0 
GO  TO  5 

3  RN=36.0 
OP=180.0 
GO  TO  5 

*»  RN=1.2566^ 
C=1.0 
D=0.0 

5  CALL  RERUN 
RETURN 

6  CALL  ENDRa(NPLOT) 
END 
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